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Abstract. Wood density is considered a key plant trait, affecting meclianical and physiological performance, yet its 
biological meaning is still rather unclear. Accordingly we investigated the anatomical underpinnings of wood density in 
trees and shrubs. We measured wood density and anatomical traits in distal stems 4-10 mm diameter under bark in 
24 Australian species. Proportions of wood components that are functionally distinct were analysed, including fibre 
wall and lumen, vessel wall and lumen, and axial and ray parenchyma. Wood density was mainly driven by the 
density of wood outside vessel lumens (densityNv) rather than by vessel lumen fraction. In turn, densityNv variation 
was chiefly affected by fibre wall and lumen fractions. Considerable anatomical variation was observed at a given 
densityNv, especially among medium-densityNv species (0.60-0.85 g cm"^); this range of medium densityNv roughly 
translates to 0.50-0.75 g cm"^ of overall density. The anatomy of these species formed a continuum from low fibre 
lumen and medium parenchyma fractions to medium fibre lumen and low parenchyma fractions. Our data suggest 
that wood density is an emergent property influenced by a complex anatomy rather than an unambiguous functional 
trait, particularly in medium-density species. With much anatomical variation, they likely represent a wide range of eco- 
logical strategies. 

Keywords: Ecological strategies; fibres; parenchyma; rays; tissue fraction/proportion/percentage/volume; vessels; 
wood anatomy. 



Introduction 

Plants vary significantly in theirecological, physiological and 
mechanical properties or 'traits' both across climate and 
even within a site (Westoby et al. 2002; Wright et al. 2004; 
Chave eta/. 2009).This indicates that there are multiple solu- 
tions to the problem of how to be a successful plant and 
maintain species continuity. Different solutions can be 
called plant ecological strategies (Westoby et al. 2002). 

Wood density has been suggested as a key player in 
plant ecological strategies (Chave et al. 2009). Firstly, 
wood density has been linked with hydraulic strategies. 



Denser woods tend to operate at more negative water 
potentials (Ackerly 2004; Bucci et al. 2004; Santiago et al. 
2004; Jacobsen et al. 2007b, 2008; Gotsch et al. 2010) 
and to have greater cavitation resistance than low-density 
woods (Hacke et al. 2001; Jacobsen et al. 2005; Pratt et al. 
2007; Lens et al. 2011). Wood density has been studied in 
relation to hydraulic conductivity but the results are incon- 
clusive, showing either a negative relationship between 
the two traits (measured conductivity, Stratton et al. 
2000; Bucci et al. 2004; Santiago et al. 2004) or no relation- 
ship (theoretical conductivity, Poorter etol. 2010; Fan etal. 
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2012). Also, density has been found to correlate negatively 
with capacitance (Meinzer et al. 2003, 2008; Pratt et ai 
2007; Scholz et al. 2007). Secondly, wood density has 
been associated with plant mechanical strategies where 
denser woods tend to be stiffer and more resistant to 
breakage at a given wood diameter (Chave et al. 2009). 
However, it has also been suggested that plants can 
build thicker stems to compensate for their lower density 
(Anten and Schieving 2010; Larjovaaro and Muller-Londau 
2010; Butler et al. 2011). Thirdly, denser woods might be 
more resistant to pathogen attacks (Augspurger and 
Kelly 1984; Romero and Bolker 2008). Wood density has 
also been widely discussed in relation to life-history strat- 
egies. For example, species with denser wood tend to ex- 
perience lower stem mortality rates (Putz et al. 1983; 
Kraft et al. 2010; Poorter et al. 2010). Growth rate is 
another important component of life history. Growth 
rate can be expected to negatively relate to wood 
density on the basis that higher investment in mass per 
volume should slow down growth (Enquist et ai 1999) 
and, while generally true, the correlation is not always 
strong (Poorter et al. 2008, 2010; Chave et al. 2009; 
Wright et al. 2010; Fan et ai 2012). Across species, wood 
density can vary with environmental factors such as tem- 
perature (Wiemann and Williamson 2002; Swenson and 
Enquist 2007; Martinez-Cabrera et ai 2009) and precipita- 
tion (Barajas-Morales 1985; Wiemann and Williamson 
2002; Swenson and Enquist 2007; Martinez-Cabrera et ai 
2009; Zhang et ai 2011), although not in all studies (Ter 
Steege and Hammond 2001; Wiemann and Williamson 
2002; Muller-Landau 2004). Despite this broad climate- 
related patterning, wood density also tends to vary quite 
widely among co-occurring species (Wiemann and 
Williamson 2002; Muller-Landau 2004). 

Thus there ore many potential functional roles for wood 
density but also a number of unresolved questions about 
each potential role. Larjavaara and Muller-Landau (2010) 
argued that some observed correlations may not be 
causal, but rather may reflect correlated selection on 
other traits. Further, if wood density has multiple func- 
tions, then it might not be a very good predictor of any 
one of them. In any event, a useful step towards resolving 
this problem is to ask what are the structural underpin- 
nings of wood density variation. The premise of the work 
reported here was that if structural underpinnings of 
wood density variation were rigorously quantified and 
better understood, this might help to explain complexities 
in the functional implications of wood density. 

Most studies of anatomical components of angiosperm 
wood density span only one or a few species (Schulz 1957; 
Taylor 1971; Taylor and Wooten 1973; Ezell 1979; Vurdu 
and Bensend 1980; Fukuzawo 1984; Bosman et ai 1994; 
Stokke and Manwiller 1994; McDonald et al. 1995; Lei 



et al. 1996; Denne and Hale 1999; Rona et ai 2009) or 
focus on commercial woods (French 1923 as cited in 
Panshin and de Zeeuw 1980; Manwiller 1973 as cited in 
Koch 1985). Fewer studies make comparisons across a 
brood number of species (Fujiwaro et ai 1991; Jacobsen 
et ai 2007a; Martinez-Cabrera et al. 2009; Poorter et ai 
2010; Fichtler and Worbes 2012). 

Wood is a complex tissue composed of three main cell 
types: vessels that transport water, fibres responsible for 
mechanical strength, and parenchyma that stores and 
transports nutrients. These tissues have different struc- 
tural characteristics and their relative proportions within 
wood influence wood density. Vessel lumens have essen- 
tially zero density; fibre and vessel walls and parenchyma 
have positive density. Vessel fraction has variously shown 
either negative or no correlation with wood density 
(Preston et al. 2006; Jacobsen et ai 2007a; Mitchell et al. 
2008; Martinez-Cabrera et ai 2009; Poorter et ai 2010; 
Zonne et ai 2010; Gleason et ai 2012). Parenchyma is 
another commonly occurring tissue, which has been 
reported to have positive, negative or no relationship 
with density (Taylor 1969; Fujiwora 1992; Jacobsen et ai 
2007o; Martinez-Cabrera et ai 2009; Rana et al. 2009; 
Poorter et ai 2010). Wood density is generally well corre- 
lated with fibre properties, especially fibre wall fraction 
(Fujiwaro et ai 1991; Jacobsen et ai 2007a; Martinez- 
Cabrera et al. 2009). However, it is unclear how these dif- 
ferent fibre traits are interrelated with each other and 
consequently how these interrelations influence wood 
density. 

Most previous work linking anatomy with density has 
concentrated on vessels with relatively little attention 
given to the other tissues. Furthermore, among the 
studies investigating all the major tissues (vessels, paren- 
chyma and fibres) only one focused on the wood of twigs in 
17 species studied by Jacobsen etal. (2007a). Twigs are im- 
portant, being in direct spatial and functional contact with 
leaves and having been commonly subjected to physio- 
logical and ecological measurements. In this paper, we 
investigate wood from twigs from a wide range of angio- 
sperm tree and shrub species growing in various environ- 
ments (24 species from four sites in eastern Australia). 
We address two main unresolved issues: (i) Which fibre 
properties have the most decisive effect on wood density 
and how ore those properties interrelated with each 
other? (ii) How do vessel and parenchyma proportions 
influence wood density? 

Methods 

Plant material and sites 

Four sites were chosen that spanned a wide range of tem- 
perature and aridity in eastern Australia [see Supporting 
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Information]. The objective of site selection was to gener- 
ate a broad range of trait values rather than to enable site 
comparisons. All carried natural, undisturbed vegetation 
growing on oligotrophic soils on flat to slightly sloping 
terrain. Two locations in Tasmania, at ~43°S, represented 
low mean annual temperature (MAT; 10.6 "C), and two 
locations in Queensland near 18°S represented higher 
MAT (c. 22.5 °C). Within each latitude, two locations were 
chosen so as to differ markedly in aridity index (AI; 
Willmott and Feddema 1992), the ratio of mean annual 
precipitation (MAP) to potential evapotranspiration (PET). 
In both Queensland and Tasmania the wetter site had an 
AI c. 1.0 and the drier site c. 0.6. MAP, MAT and PET were 
obtained from GIS (geographic information system) 
layers from the Australian Bureau of Meteorology. 

At each of the four sites, six abundant and phylogenet- 
ically distinct woody eudicot species were chosen for sam- 
pling (species listed in Table 1). One species was sampled 
at two sites, yielding a total of 23 species from eight fam- 
ilies. Distal, sun-exposed twigs of trees and shrubs were 
collected from three replicate individuals per species. 
The diameter under bark of twigs varied from 4 mm in 
plants with little pith to 10 mm in plants with higher pith 
content. Consequently the diameter of wood, excluding 
bark and pith, was 4-5 mm. This plant material is referred 
to here as twigs, although in several small shrub species 
'twigs' were the main stems. Plant material was cut into 
segments 10-15 cm long and kept wet in sealed plastic 
bags in the refrigerator (4 °C). Wood density was measured 
within a week from collection; other parts of the same 
twigs were placed in fixative for later measurement of 
anatomical properties (details below). 

Wood density 

Wood density was measured on segments 3-5 cm long 
for each twig sample. Bark and pith were removed and 
measurements were carried out on xylem only. In this 
paper we refer to xylem as 'wood'. After removing bark 
and pith, wood pieces were soaked in water for at least 
48 h prior to volume measurement. Then a beaker filled 
with water was placed on a balance (0.0001 g, Mettler 
AE 160). A thin wire platform was suspended in water so 
that it did not touch the side or bottom of the beaker. 
The balance was tared before each measurement and a 
sample was gently placed on the platform. The mass of 
displaced water was read from the balance. From stand- 
ard water density of 1 g cm"^ and knowing the mass of 
displaced water, we calculated sample volumes applying 
Archimedes' buoyancy principle (e.g. 1 g of displaced 
water equals 1 cm^ volume). Samples were then placed 
in paper envelopes and dried at 70 °C for at least 72 h. 
Wood density was calculated as the dry mass divided by 
water saturated volume (g cm"^). 



Table 1. Sampled sites, species names and families. 



Site 


Species name 


Family 


Cool-wet 


Allocasuarina monilifera 


Casuarinoceae 




Aotus ericoides 


Fabaceae 




Banksia marginata 


Proteaceae 




Eucalyptus amygdaima 


Myrtaceae 




Leptospermum scoparium 


Myrtaceae 




Leucopogon ericoides 


Ericaceae 


Cool-dry 


Bossiaea cinerea 


Fabaceae 




Davesia latifolia 


Fabaceae 




Epacris impressa 


Ericaceae 




Eucalyptus tenuiramis 


IMyrtoceae 




Leucopogon ericoides 


Ericaceae 




Persoonia juniperina 


Proteaceae 


Hot-wet 


Acacia mangium 


Fabaceae 




Allocasuarina torulosa 


Casuarinoceae 




Alphitonia excelsa 


Rhamnoceae 




Chionanthus ramiflorus 


Oleoceae 




Eucalyptus platyphylla 


IMyrtoceoe 




Ixora timorensis 


Rubioceoe 


Hot-dry 


Acacia flavescens 


Fabaceae 




Corymbia intermedia 


IMyrtoceoe 




Gastrolobium grandiflorum 


Fabaceae 




Greviilea parallela 


Proteaceae 




Lopliostemon suaveolens 


IMyrtoceoe 




Persoonia falcata 


Proteaceae 



Anatomy 

To obtain anatomical cross-sections the material was first 
fixed in formalin-acetic acid-alcohol (FAA) for 4 weeks. 
The FAA was prepared in proportions of 5 : 5 : 90 (forma- 
lin: glacial acetic acid: 70% ethonol; Gerlach 1972). 
After 4 weeks, the fixative was replaced with 70 % 
ethanol. This 70 % ethonol was then replaced two more 
times within 10 days to further wash the fixative out. 
Final replacement of alcohol was used as a long-term 
storage medium. Segments for image analysis were rehy- 
drated by immersion in 50 % ethanol and after 2-5 days 
in 30 % ethanol. Cross-sections were cut with a sledge 
microtome (Reichert, Vienna, Austria) at 10-20 |xm thick- 
ness using disposable blades (Model A35, Feather Safety 
Razor Co. Ltd, Japan). For better contrast and tissue identi- 
fication, sections were stained with safronin 0 (Gurr 
Microscopy Materials, BDH Chemicals Ltd, UK) for lignified 



AoB PLANTS www.aobplants.oxfordjournals.org 



©The Authors 2013 



3 



Ziemihska et at. — Anatomical underpinnings of wood density 



cell walls (10 mln) and with Janus green B (Gurr's, London, 
UK) for cytoplasm (10 min). Sofranin 0 solution used 2 g of 
stain in 100 mL of distilled water (Ruzin 1999) and Janus 
green B used 0.1 g of stain and 1 mL of glacial acetic acid 
in 100 mL of distilled water (Conn et al. 1960). Sections 
were rinsed in distilled water after each staining session. 
Afterwards, they were mounted in glycerol on a slide, 
covered with a cover slip and sealed with noil polish. Mea- 
surements were mode only on cross-sections, but to assist 
in interpreting and identifying cell types, longitudinal tan- 
gential and radial sections and macerations were also 
made. For macerations, small shavings were placed in 
vials filled with Franklin's solution: glacial acetic acid and 
6 % hydrogen peroxide in proportions of 1 : 1 (Franklin 
1945). These vials were loosely covered with Porafilm 
tope and heated in the oven at 60 °C for 1-2 days. 
Tissues were then rinsed with distilled water, stained 
with sofranin 0 (10 min) and gently squashed onto a 
microscope slide. 

Microphotogrophsofcross-sections were token at xlOO 
and x400 magnifications using a digital camera (Scion 
Corporation, CFW-1310C, USA) attached to a light micro- 
scope (Olympus BX 50F, Olympus Co. Ltd, Japan) and 
image capturing software Scion Visicopture, version 1.4 
(Scion Corporation). Two to three images of the same 
area at different focal planes were token and stacked in 
Photoshop CS4 (Adobe Systems Incorporated, USA). Cross- 
sections were bigger than the field of view; therefore, 
dozens (for xlOO) or around 10 (for x400) images per 
cross-section were taken and then merged in Photoshop, 
giving rise to images of whole cross-sections at xlOO 
and of one narrow transect at x400. Tissue cross- 
sectional areas and vessel traits were measured on one 



wedge-shaped transect per replicate (xlOO; Fig. 1) and 
fibre characteristics (x400) on one rectangular transect, 
both stretching from pith to cambium. The radial transects 
were chosen to be the most representative for a section 
and the tension wood was avoided where possible. The 
transect borders were approximately parallel to the rays 
and followed middle lamella so that no open cells were 
positioned on the borders. Tissue areas, vessels and fibre 
walls were manually coloured in Photoshop using a 
Cintiq 21UX graphic tablet (Wacom Co., Ltd, Japan). Proto- 
xylem and newly produced xylem were excluded from 
analysis. Larger regions were inspected in species with 
larger vessels or with more variable structure in the tan- 
gential direction. The measured area ranged from 0.21 
to 0.69 mm^ for species from cool sites and from 0.31 to 
1.9 mm^ for species from hot sites. On average, for each 
image there were 170 vessels measured ( + 145) of all 
sizes including vessel tails. Fibre, fibre wall and lumen 
areas were measured for an overage of 170 fibres per 
sample ( + 57) lying in two parallel rows from pith to 
cambium. Colour-coded images were analysed with 
Image-Pro Plus version 2.0.0.260 (Media Cybernetics, 
Inc., USA; Fig. 1). 

Proportions of all major wood cell types (vessel wall and 
lumen, fibre wall and lumen, axial parenchyma, rays and 
trocheids) and their properties (mean vessel lumen area 
and mean fibre lumen and wall area) were quantified. 
Here, mean vessel lumen area is also called 'mean vessel 
area' for brevity. Mean vessel area, fibre lumen area and 
wall area were calculated as the arithmetic mean across 
all relevant cells (vessels or fibres) measured within a 
given radial sector of a sample. In addition to these 
traits we calculated theoretical maximum hydraulic 




Figure 1. A twig cross-section of Grevitlea paraileta, Proteaceae. The radial sector of the stained section is shown on the left side and the processed 
image on the right. Colours in the processed image denote different tissue types: blue, vessel lumen; purple, vessel wall; green, rays; orange, axial 
parenchyma; brown, fibres. The scale bar corresponds to 100 |j.m. 
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conductivity (called here 'theoretical conductivity'; also 
known as 'potential conductivity'). Conductivity for the 
lumen of each conduit was calculated from the Hagen- 
Poiseuille expression (-Trr^j/CSix), where r represents 
mean vessel lumen radius and |x represents water 
dynamic viscosity (assuming a standard viscosity value 
of 1 at 20 "C). The sum of all vessel lumen conductivities 
per unit cross-section is theoretical conductivity. This 
quantity expresses conductivity variation across species 
due to conduit diameters, but it should be thought of as 
the theoretical maximum since it does not capture 
effects of end-wall resistances, blockage by tyloses, 
effects of emboli, etc. 

Traditionally, wood is considered a complex tissue com- 
posed of several cell types (Evert 2006). However, we refer 
to those cell types as 'tissues' for brevity and also because 
they perform distinctly different functions. Anatomical ter- 
minology follows the 'lAWA list of microscopic features for 
hardwood identification' (lAWA Committee 1989). Vascular 
and/or vosicentric tracheids occurred in 13 species and are 
referred to hereafter as tracheids. Tracheids were first 
determined in macerated wood and then identified on a 
cross-section on the basis of the number of pits and size 
of the pit border (both resembling that of vessel pits) and 
cell size (lAWA Committee 1989; Sano et al. 2011). Axial 
apotracheol and paratracheal parenchyma were collective- 
ly measured as axial parenchyma. None of the species 
studied here had storied rays. Tissue types were expressed 
as the fraction of a tissue per cross-sectional area (Fig. 1). 
Tissue fractions of the area outside vessel lumen (non- 
vessel area) were calculated as tissue fraction multiplied 
by non-vessel area fraction. The term 'non-vessel' area is 
used for brevity and it includes vessel walls, fibre walls 
and lumens, axial and ray parenchyma, and tracheids. 
Non-vessel-lumen quantities ore hereafter denoted by 
subscript 'NV, e.g. fibre wall fractionNv, wood densityNv, 
etc. Fibre wall proportion in a given fibre was expressed 
as a proportion of the total fibre area. Fibre wall fraction 
was obtained by multiplying mean fibre wall proportion in 
a fibre by the fibre fraction per cross-section. Fibre lumen 
fraction was similarly calculated from mean fibre lumen 
proportion in a fibre multiplied by fibre fraction. 

Statistical analysis 

We collected measurements for 23 species, one of which 
occurred at two sites and was considered as two entities, 
giving a total of 24 data points analysed. Measurements 
were carried out on three replicate individuals per 
species and the trait values were averaged for compari- 
sons across species. Wood density, vessel lumen, sum of 
ray and axial parenchyma, ray parenchyma and fibre 
wall fractions of total wood area as well as non-vessel 
wood area were all approximately normally distributed 



(Shapiro- Wilk test, P<0.05). Vessel wall, axial paren- 
chyma, tracheids and fibre lumen fractions of total wood 
area as well as of non-vessel wood area were right-skewed 
and were transformed to generate approximately normally 
distributed variables. Logio transformations normalized 
all distributions with the exception of tracheid fraction 
and tracheid non-vessel fraction. We used ordinary 
least-squares regression to assess bivariate relationships 
(SigmaPlot, Systat, San Jose, CA, USA). 

Results 

Wood density and tissue proportions varied significantly 
across species as illustrated in Fig. 2 [see Supporting Infor- 
mation]. Wood density varied more than 2-fold, from 0.37 
to 0.83 g cm"^. Figure 2 shows tissue fractions averaged 
across all species (bar at the top) and for each individual 
species separately (the remaining bars). The mean fibre 
fraction was 0.52 + 0.09 (hereafter numbers represent 
overage fraction + one standard deviation). Fibre varied 
approximately 2-fold across species and was the most abun- 
dant tissue type. Fibre fraction could be partitioned into fibre 
walls (0.45 + 0.08, 2-fold variation; brown bars in Fig. 2) and 
fibre lumens (0.08 + 0.07, c. 60-fold variation; yellow bars). 
On average, parenchyma occupied 0.25 of wood cross- 
sectional area and varied almost 3-fold across species. 
It consisted of axial parenchyma (0.10 + 0.05, 6-fold 
variation) and ray parenchyma (0.15 + 0.05, 4.5-fold 
variation). Vessels occupied 0.20, where 0.15 + 0.03 con- 
sisted of lumens (varying c. 2-fold) and 0.05 + 0.03 of 
vessel walls (varying 4.5-fold). Tracheids occurred in 
just 13 of 24 species and occupied only small fractions 
(0.02 + 0.03 averaged across all 24 species, 8.5-fold 
variation across the 13 species that had tracheids). 

Alternatively, wood can simply be divided into two com- 
ponents: vessel lumen fraction and non-vessel fraction, 
which encompasses all tissues other than vessel lumens. 
The density of the non-vessel fraction and tissue fractions 
within the non-vessel fraction ore indicated hereafter by 
the subscript 'NV, e.g. densitynv, fibre fractionNv- Since 
vessel lumen has zero density, overall wood density is (by 
definition) the product of the non-vessel fraction density 
(densitynv) and the non-vessel fraction itself (fractionNv): 
density = densitynv x fractionNv (Preston et al. 2006; Zonne 
et al. 2010). These three quantities con be log transformed 
and the equation then becomes a sum: log(density) = 
log(densityNv) -I- log(fractionNv)- Hence, when densityNv 
is plotted against fractionNv on log -log axes (Fig. 3), iso- 
lines of resulting overall wood density can be constructed 
and used to aid interpretation. Variation across species in 
densityNv (y-oxis in Fig. 3) was four times greater than vari- 
ation in fractionNv (x-oxis), and thus in this species set was 
0 for stronger determinant of variation in overall wood 
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Tissue fractions 
0.4 0.6 




I Fibre wall 
I Vessel wall 
I Tracheids 



I I Axial parenchyma 
1 I Vessel lumen 
I I Fibre lumen 



Average 

Bossiaea cinerea, Fabaceae 
Acacia flavescens, Fabaceae 
Eucalyptus tenuiramis, Myrtaceae 
Leptospermum scoparium, Myrtaceae 
Leucopogon ericoides', Ericaceae 
Gastrolobium grandiflorum, Fabaceae 
Epacris impressa, Ericaceae 
Aotus ericoides, Fabaceae 
Leucopogon ericoides, Ericaceae 
Persoonia jur)iperina, Proteaceae 
Daviesia latifolia, Fabaceae 
Corymbia intermedia, Myrtaceae 
Persoonia falcate, Proteaceae 
Grevillea parallela, Proteaceae 
Allocasuarina torulosa, Casuarinaceae 
Eucalyptus amygdalina, Myrtaceae 
Allocasuarina monilifera, Casuarinaceae 
Chionanthus ramiflorus, Oleaceae 
Lophostemon suaveolens, Myrtaceae 
Ixora timorensis, Rubiaceae 
Banksia marginata, Proteaceae 
Eucalyptus platyphylla, Myrtaceae 
Acacia mangium, Fabaceae 
Alphitonia excelsa, Rhamnaceae 



Ray parenchyma 

Figure 2. Tissue fractions for 24 species arranged in order of increasing wood density (from bottom to top). Large wood density numbers indicate 
total wood density whereas small numbers indicate non-vessel density (density mv). Mean tissue fractions across all species are shown in the bar at 
the top of the figure. Leucopogon ericoides occurred in two sites and is treated here as two separate entities. An asterisl< denotes L. ericoides from 
the cold-wet site as opposed to L. ericoides from the cold-dry site with no indicator. 



Fraction^y 
0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 




Figure 3. Relationship between fraction of wood outside vessel 
lumens (fractionfjv) and the density of that non-vessel fraction 
(wood densityNv) among 24 Australian species. FractionNv = 1 - 
vessel fraction. Each circle represents a different species (mean 
value from three replicates). Diagonal isolines represent contours of 
overall wood density, which increases towards the upper right. All 
axes are log scaled. 

density (direction across the isolines). Not surprisingly 
then, densityNvond overall density were tightly correlated 
with each other (r^ = 0.95, P < 0.001). Vessel lumen frac- 
tion was only loosely (negatively) correlated with overall 



wood density (r^ = 0.20, P = 0.027). Therefore, the follow- 
ing analyses concentrate entirely on densityNvond its ana- 
tomical components [see Supporting Information]. 

Total fibre fractionNv (fibre lumen fractionNv plus fibre wall 
fractionsNv) was unrelated to densityNv (Fig. 4). Species with 
the some fibre fractionNv varied widely in wall proportion 
relative to lumen within a fibre (as indicated by the width 
of the 'donut' rings in Fig. 4). Figure 4 shows that lowest- 
density (<0.5 g cm"^) species had high fibre fractionNv 
but their fibres had low wall proportion (lower right of the 
graph). High-densityNv species (>0.85 g cm"^) also hod 
high fibre fractionNv but their fibres hod large wall proportion 
(upper right of the graph). A substantial number of species 
located in the middle of the graph with medium densityNv 
(0.60-0.85 g cm"^) had variable fibre fractionNv arid 
fibre wall proportion within a fibre. Also, fibre wall proportion 
within a fibre was positively correlated with wood densityNv 
(r^ = 0.62, P < 0.001). 

DensityNv was positively correlated with fibre wall 
fractionNv (f^ = 0-^0; Fig. 5A), and negatively with fibre 
lumen fractionNv (("^ = 0-56; Fig. 5B). Since the majority 
of species had low fibre lumen fractionNv (i-e. the data 
were right-skewed), the two species with highest fibre 
lumen fractionNv had a strong influence on these relation- 
ships (lower left of Fig. 5A, lower right of Fig. 5B). That said, 
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the correlations were still present across the other species 
considered on their own (r^ = 0.26 and P = 0.014, = 
0.18 and P = 0.049, respectively). The species with the 
lowest fibre wall fractionNv (upper left in Fig.SA) wasDov/e- 
sia latifoUa, which had a high amount (fraction of 0.07) of 
thick-wall trocheids. Presumably, this high fraction of tra- 
cheid wall contributed to the quite high densityNv of this 
species (high, given its very low fibre fraction). 

Next we asked how fibre wall and lumen fractionsNv 
were related to one another. Figure 6 shows on approxi- 
mately triangular relationship between fibre wall fractionNv 




0.5 0.6 0.7 0.8 

Fibre fractionNv 

Figure 4. Relationship between non-vessel density (wood densityNv) 
and fibre fraction in non-vessel fraction (fibre fractionfjv). Each'donut' 
circle symbolizes one species. The width of the donut border (block) 
represents fibre wall proportion within on individual fibre and the 
width of the hole (white) represents fibre lumen proportion. These 
proportions were estimated from individual fibres (as fibre wall 
area — or lumen area — divided by total fibre area), for 75-314 fibres 
per replicate (mean 170), and then across three replicates per species. 



and fibre lumen fractionNv- The highest-densityNv species 
(>0.85gcm"^; four large symbols, upper left of Fig. 6) 
had high fibre wall fractionNv and low fibre lumen 
fractionNv- Medium-densityNv species (0.60-0.85 g cm"^; 
18 medium-sized symbols in Fig. 6) had variable fibre wall 
fractionNv ond variable fibre lumen fractionNv The two 
lowest-densityNv (<0.50 g cm"^) species had the highest 
fibre lumen fractionNv and low fibre wall fractionNv (two 
smallest symbols, lower right of Fig. 6). The categories of 
high-, medium- and low-density species are used here for 
easy reference, but in fact, the trait values are continuous 
and no clear boundaries can be indicated. The species 
with lowest fibre wall fractionNv (lower left in Fig. 6) was 
D. latifolia (see the comment about trocheids above). 

The second most abundant tissue, parenchyma, was 
not correlated with densityNv nor were its components, 
rays and axial parenchyma (all P > 0.7). Similarly, vessel 
wall fractionNv was unrelated to densityNv Nevertheless, 
both parenchyma and vessel wall fractionsNv indirectly 
affected densityNv Figure 6 illustrates that there was con- 
siderable variation in densityNv (indicated by symbol size) 
at a given fibre wall fractionNv especially at lower wall 
fractionNv (lower half of the graph), and considerable vari- 
ation in densityNv at a given fibre lumen fractionNv espe- 
cially at lower lumen fractionNv (left half of the graph). 
These variations in densityNv could be partially explained 
by parenchyma and vessel wall fractionsNv At a given 
fibre wall fractionNv densityNv was positively correlated 
with parenchyma and vessel wall fractionsNv (r^ = 0.15, 
P = 0.064 and = 0.29, P = 0.007, respectively) and nega- 
tively with fibre lumen fractionNv (f"^ = 0.45, P < 0.001; all 
relationships tested against residuals from a regression 
of wood densityNv on fibre wall fractionNv)- Conversely 
at a given fibre lumen fractionNv (i-S- tested against 
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Figure 5. Relationships between non-vessel density (wood densityNv) and (A) fibre wall fraction in non-vessel fraction (fibre wall fractionNv) and (B) 
fibre lumen fraction in non-vessel fraction (fibre lumen fractionNv). Each circle represents a different species (mean value from three replicates). 
***P< 0.001. 



AoB PLANTS www.aobplants.oxfordjournals.org 



©The Authors 2013 



7 



Ziemihska et at. — Anatomical underpinnings of wood density 



residuals fronn a regression of wood densityNv on fibre 
lumen fractionNv), wood densityNv was negatively corre- 
lated with parenchyma fractionfgy (r^ = 0.28, P = 0.008), 
not correlated with vessel wall fractionNv (P = 0.29) and 
positively correlated with fibre wall fractionNv ir^ = 0.40, 
P< 0.001). Additionally, parenchyma fracti onivjv was 
tightly negatively correlated with total fibre fractionNv 
(r^ = 0.74, P< 0.001). The only remaining tissue, trach- 
eids, occupied on average a very small fractionNv and 
was not subjected to detailed analysis. 

Mean vessel area and theoretical conductivity [see 
Supporting Information] were negatively correlated 
with overall wood density across the species studied (r^ = 
0.225, P = 0.019 and = 0.33, P= 0.0036, respectively, 
Fig. 7A and B). However, these relationships were mainly 
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Figure 6. Relationship between the fraction of the non-vessel frac- 
tion that is fibre wall or is fibre lumen. Each symbol represents a differ- 
ent species (mean value from three replicates). Symbol diameter is 
proportional to non-vessel wood density (wood densityNv), with the 
biggest symbol indicating highest densityNv (0.92 g cm~^) and the 
smallest symbol indicating lowest densityNv (0.44 g cm"^). 



driven by a few species with particularly large vessel 
lumens and theoretical conductivity (Fig. 7A and B). Among 
the majority of species there was considerable variation in 
mean vessel area and conductivity at a given wood density, 
and little relationship between the two (Fig. 7A and B). 

Discussion 

This study aimed to describe the anatomical components 
of wood density in twigs of 24 Australian tree and shrub 
species. Properties of fibres, the most abundant tissue, 
had the strongest effect on wood density variation, as 
has been shown in some previous studies (Fujiwara et al. 
1991; Jacobsen et ai 2007a; Martinez-Cabrera et ai 
2009). However, our results contribute to a more compre- 
hensive understanding of wood structure and its influence 
on density in twigs. Here, we discuss in detail the proper- 
ties of fibres and other tissues as components of density. 

Wood density and its anatomical components 

Density of tissue outside vessel lumens (densityNv), rather 
than vessel lumen fraction, was the main driver of overall 
wood density variation. This result agrees with a compari- 
son made across 584 species that considered main stem 
wood (Zanne et al. 2010), implying that densityNv deter- 
mines overall density in both twigs (this study) and main 
stems similarly. Accordingly, the discussion here is direc- 
ted towards densityNv ond individual tissue fractions 
within the non-vessel part of the wood (indicated by sub- 
script 'NV, e.g. fibre wall fractionNv)- We compare these 
results with the results for overall wood density reported 
by other studies, on the basis that overall wood density 
and densityNv are closely correlated (r^ = 0.95, P < 
0.001, this study). 

All tissue fractionsNv influence densityNv but fibre 
wall and lumen fractionsNv are the most important. The 
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strong influence of fibre wall fractionNv was due both to its 
high proportion and its variability, while the influence of 
fibre lumen fractionNv was associated chiefly with its 
high variability (c. 60-fold). Other studies have consistently 
found a positive relationship between fibre wall fraction 
and density: in trunk wood among 50 Japanese trees 
(Fujiwara et al. 1991) and 61 North and South American 
shrubs (Martinez-Cabrera et al. 2009); and in twig wood 
of 17 South African shrubs (Jocobsen et al. 2007a). Fibre 
lumen fraction has received less attention but has also 
been shown (in concordance with our study) to have a 
negative relationship with wood density (Martinez- 
Cabrera et al. 2009; Rana et al. 2009). The second most 
abundant tissue, parenchyma, did not correlate with 
densityNv in our study nor in tree and shrub trunks 
(Martinez-Cabrera et al. 2009; Poorter et al. 2010; Fichtler 
and Worbes 2012) but correlated negatively with twig 
wood density across 17 species (Jacobsen et al. 2007o). 
These discrepancies might be caused by variable densities 
of parenchyma tissue itself (Taylor 1969; Fujiwara 1992; 
Guilley and Nepveu 2003) or by varied correlations 
between parenchyma and fibre wall and lumen fractions 
(see below). Plausibly, these discrepancies could also 
stem from different relationships between ray and axial 
parenchyma. We did not find any cross-correlation 
between those two components of parenchyma nor 
were they individually related to density. In contrast, 
Martinez-Cabrera et al. (2009) reported that ray and axial 
parenchyma were negatively correlated with each other 
and individually correlated with density (axial paren- 
chyma positively correlated, rays — negatively). In that 
study these links were strongly associated with environ- 
mental variables (MAP, MAT and Al). These findings imply 
that different functional trade-offs can be related to ray 
and axial parenchyma individually and also that their 
link is affected by climate. Possibly these trade-offs may 
be more pronounced in trunk wood, as opposed to the 
twig wood studied here. 

We did not find a direct relationship between paren- 
chyma and densityNv. Nevertheless, our results imply 
that parenchyma together with fibre lumens con influence 
densityNv variation in a less direct way. At a given fibre wall 
fractionNv, densityNv depended on the parenchyma 
fractionNv relative to fibre lumen fractionNv- DensityNv 
was marginally positively correlated with parenchyma 
fractionNv and negatively with fibre lumen fractionNv Par- 
enchyma has higher density than fibre lumen, which has 
zero density. Therefore, higher parenchyma fractionNv 
relative to fibre lumen fractionNv increases densityNv, 
and conversely higher fibre lumen fractionNv decreases 
densityNv We note that parenchyma fractionNv was only 
weakly correlated with densityNv at a given fibre wall 
fractionNv This weak correlation possibly stems from 



variable parenchyma tissue densities (Taylor 1969; 
Fujiwara 1992; Guilley and Nepveu 2003). 

Fibre fractionNv (sum of fibre wall and lumen fractionsNv) 
was not associated with densityNv because of wide vari- 
ation in wall proportion within a fibre (top right to low 
right in Fig. 4). Poorter et al. (2010) suggested a similar ex- 
planation, but as for as we are aware, this issue has not 
been quantitatively clarified until now. 

Anatomical traits that are not expressed as fractions of 
wood volume have less direct effects on density. We found 
that densityNv was correlated with fibre wall proportion 
within a fibre. Previous studies hove also shown that 
density can be related to fibre lumen diameter (Jacobsen 
etal. 2007o; Martinez-Cabrera eto/. 2009), fibre wall thick- 
ness (Fujiwara et al. 1991; Martinez-Cabrera et al. 2009) 
and fibre wall to lumen ratio (Martinez-Cabrera et al. 
2009). We believe those traits would deliver a more in- 
sightful understanding of wood density when combined 
with information about fractions, e.g. the relationship 
among densityNv. fibre fractionNvand fibre wall proportion 
in a fibre, as described above. Similarly for mean vessel 
area, which in this study was weakly negatively correlated 
with wood density. Mean vessel area per se should not 
affect wood density. Rather, vessel fraction (vessel area 
multiplied by vessel number per area) should be causally 
linked with density variation. 

Variability of anatomical structures 

The discussion so far has focused on wood density variation. 
However, we also found considerable anatomical variation 
within a given range of density. Hereafter, we use the arbi- 
trary categories of 'medium', 'high' and 'low' density and 
tissue fractions. However, we observed a continuum of 
trait values and the categories are used only for conveni- 
ence. Species with medium densityNv (0.60-0.85 g cm"^) 
showed broader structural variability than high- and 
low-densityNv species (>0.85 and <0.5 g cm"^, respect- 
ively). The concept is illustrated in Fig. 8 and examples 
of cross-sections approximately corresponding to Fig. 8 
ore shown in Fig. 9. Species with high densityNv (large 
symbols in Fig. 6, top corner in Figs 8 and 9A) had 
the highest fractionNv of fibre wall and small fibre 
lumen fractionNv Their total fibre fractionNv was high 
and parenchyma fractionNv was low. In contrast, medium- 
densityNv species had more variable fibre wall, fibre lumen 
and parenchyma fractionsNv (medium symbols in Fig. 6, 
middle in Figs 8 and 9B and C). Consequently, a spectrum 
of possible architectures may be outlined where the two 
ends of the spectrum are (i) low fibre, low fibre lumen 
and high parenchyma fractionsNv (middle left in 
Figs 8 and 9B) and (ii) high fibre, medium fibre lumen 
and low parenchyma fractionsNv (middle right in Figs 8 
and 9C). The lowest-density species (<0.5 g cm"^, small 
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Figure 8. A scliematic representation of tlie grapin in l^ig. 6 flipped 
clocl<wise by 45 \ Tlie actual graph from Fig. 6 is repeated here in an 
inset and also flipped 45° clockwise. The diagram represents four 
cross-sections of potential anatomical structures in low-, medium- 
ond high-density species. Each hexagon within the three squares 
indicates a fibre cell consisting of dark fibre wall and bright fibre 
lumen. The green area on the right of each square indicates paren- 
chyma. Wood densityNv increases towards the top of the diagram. 

symbols in Fig. 6, bottom comer in Figs 8 and 9D) in this 
study had high fibre lumen fractionNv- However, it is pos- 
sible that low-density wood could also be composed of 
large parenchyma fractionSfgv and small fibre lumen 
fractionsNv More studies are needed to clarify the range 
of anatomies exhibited by low-density species. To our 
knowledge, we are the first to set out this triangular 
scheme relating wood anatomy to density. Variability of 
structures in medium-density species implies that there 
may be a wider range of ecological strategies available 
to these species. 

Ecological and physiological considerations 

It has been shown that denser woods tend to operate at 
more negative minimum water potentials (Ackerly 2004; 
Bucci et al. 2004; Santiago et al 2004; Jacobsen et al. 
2007b, 2008; Gotsch et al. 2010) and have higher 



cavitation resistance (Hocke et al. 2001; Jacobsen et al. 
2005; Pratt et al. 2007; Lens et al. 2011). However, we 
and previous literature hove shown that wood density is 
most strongly related to fibre properties, not vessel prop- 
erties. Thus, the relationship among wood density, 
minimum water potential and cavitation resistance does 
not appear to be directly influenced by wood density, but 
rather by a third, unexplained factor. Another element of 
hydraulic strategies is capacitance, which tends to be 
negatively correlated with wood density (Meinzer et al. 
2003, 2008; Pratt et al. 2007; Scholz et al. 2007). Our 
results imply that capacitance water could be stored in 
fibre lumen or in parenchyma but the exact mechanism 
of this water storage and release is not known, nor is it 
clear what might be the advantage or disadvantage of 
fibre lumen capacitance compared with parenchyma cap- 
acitance. Density also plays a role in the mechanical be- 
haviour of wood. It is a good predictor of mechanical 
strength and elasticity at a given wood diameter (Chave 
et al. 2009). However, it has also been shown that plants 
can compensate for mechanically weak wood by building 
thicker stems (Anten and Schieving 2010; Larjavaara and 
Muller-Landau 2010; Butler et al. 2011). Thick, low-density 
stems could be mechanically strong but their mainten- 
ance costs would be higher due to large surface area 
(Larjavaara and Muller-Landau 2010). Additionally, our 
results suggest that in similar density woods maintenance 
costs per sapwood cross-sectional area could vary de- 
pending on wood anatomical structure. Namely, woods 
with high fraction of living parenchyma would presumably 
incur higher maintenance costs than woods with high 
fraction of non-living fibres. What then would be the ad- 
vantage of high parenchyma fraction? Hypothetically, 
high parenchyma fraction could ensure high nutrient 
storage capacity, which may be beneficial in some eco- 
logical strategies. This issue and other questions about 
causal links between wood anatomical structure and eco- 
logical strategies require further investigation and open 
interesting pathways for future research. 

Twig and main trunk comparisons 

Caution is needed when comparing results from twigs 
with results from main trunks. It is not proven that the rela- 
tionships of wood density and anatomy are the same at 
the twig level as at the trunk level. Taylor and Wooten 
(1973) found that in five species vessel and fibre fractions 
shifted in the some direction with plant height, but this 
was not the case for ray fraction. Other studies examining 
association between vessel lumen fraction and wood 
density across a large number of species showed no rela- 
tionship in the main trunks (Martinez-Cabrera et al. 2009; 
Poorter et al. 2010; Zanne et al. 2010) but a negative rela- 
tionship in twigs (Preston et al. 2006; Jacobsen et al. 2007a; 
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Figure 9. Cross-sections tlirough twigs of four species. Tlie triangle arrangement corresponds to the one in Fig. 8 and so the cross-sections are 
examples of respective anatomies drown in that figure. (A) Bossiaea cinerea (Foboceae, wood densityMvO.92 g cm~^), (B) Grevillea parallela (Pro- 
teaceoe, wood densityNv 0.74 g cm^^), (C) Corymbia intermedia (Myrtaceae, wood densityNv 0.78 g cm^^), (D) Alphitonia excelsa (Rhamnaceae, 
wood densityMv 0.44 g cm" ^). V, vessel; FW, fibre wall; FL, fibre lumen; A, axial parenchyma; R, ray parenchyma. Fibre lumens are indistinct in (A) due 
to their small sizes. 



Mitchell et al. 2008; Gleason et al. 2012). Such disparities 
indicate that the relationship between wood density and 
tissue fractions may conceivably be different in main 
trunks than in twigs, yet it is unclear how generally this is so. 

Conclusions 

Wood density has been proposed as a key plant functional 
trait and is related to ecological strategies (Chave et al. 



2009) but relatively little is known about the underpinnings 
of these relationships. This study clarifies the anatomical 
basis for wood density variation across species. It shows 
that wood density depends on anatomical structure, but 
also that a range of very different structures can result in 
very similar wood density, especially among species with 
medium densityNv (here, 0.60-0.85 g cm"^). This conclu- 
sion suggests that there may be a wider range of ecological 
strategies among such species. Taken together, these 
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findings imply that twig wood density should not be consid- 
ered as an unambiguous indicator of plant ecological strat- 
egies. We hope this research will enhance the interpretation 
and design of ecological studies related to wood density. 
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